X-ray Linear Dichroism in cubic compounds: the case of Cr 3+ in MgAl 2 04 
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The angular dependence (x-ray linear dichroism) of the Cr K pre-edge in MgAl2 04:Cr 3+ spinel is 
measured by means of x-ray absorption near edge structure spectroscopy (XANES) and compared 
to calculations based on density functional theory (DFT) and ligand field multiplet theory (LFM). 
We also present an efficient method, based on symmetry considerations, to compute the dichroism 
of the cubic crystal starting from the dichroism of a single substitutional site. DFT shows that 
the electric dipole transitions do not contribute to the features visible in the pre-edge and provides 
a clear vision of the assignment of the ls^>3d transitions. However, DFT is unable to reproduce 
quantitatively the angular dependence of the pre-edge, which is, on the other side, well reproduced 
by LFM calculations. The most relevant factors determining the dichroism of Cr K pre-edge are 
identified as the site distortion and 3d-Sd electronic repulsion. From this combined DFT, LFM 
approach is concluded that when the pre-edge features are more intense than 4 % of the edge jump, 
pure quadrupole transitions cannot explain alone the origin of the pre-edge. Finally, the shape of 
the dichroic signal is more sensitive than the isotropic spectrum to the trigonal distortion of the 
substitutional site. This suggests the possibility to obtain quantitative information on site distortion 
from the x-ray linear dichroism by performing angular dependent measurements on single crystals. 

PACS numbers: 61.72.Bb, 78.70.Dm, 71. 15. Mb 



I. INTRODUCTION 



Transition metal elements play an essential role in 
physics (magnetic materials, superconductors...), coor- 
dination chemistry (catalysis, metalloproteins) or geo- 
physics (3d elements are major constituents of the Earth 
and planets). To understand the properties that transi- 
tion elements impart to the materials they are inserted in, 
X-ray Absorption Near Edge Structure (XANES) spec- 
troscopy has been widely used, since it provides unique 
information on their local surrounding and electronic 
structure. In particular, the position and intensity of the 
localized transitions observed at the K pre-edge (ls^3d 
transitions) are sensitive to the cation oxidation state, 
the geometry of its environment (coordination number 
and symmetry), and the degree of admixture between 
p and d orbitals. For example, the shape and area of 
the pre-edge are commonly used to quantify the redox 
states of transition elements in crystals, glasses and co- 
ordination complexes, by comparison to those recorded 
on reference compounds. 1-3 However, it is not straight- 
forward to obtain this kind of information on single crys- 
tals. Indeed, it is well known that the XANES spectra 
of non-cubic crystals show an angular dependence, when 
the polarization and the direction of the incident x-ray 
beam (here, designated as unit vectors, i and k, respec- 
tively) are varied. For cubic crystals, the problem may 
seem at first sight more simple. Electric dipole transi- 
tions (e.g., Is— >p transitions) are isotropic. They con- 
tribute mainly to the edge, but also to the pre-edge if 



one of the three following situations is encountered : (i) 
there is p-d intrasite hybridization (e.g., the crystallo- 
graphic site does not show an inversion center), (ii) the 
thermally activated vibrations remove the inversion cen- 
ter, (iii) there is p-d intersite hybridization (in samples 
highly concentrated in the investigated element). Elec- 
tric quadrupole transitions are anisotropic and the cubic 
crystal thus shows an angular dependence. The infor- 
mation carried by the pre-edge features can be derived 
for cubic crystals by taking advantage of this angular 
dependence. In particular, the respective proportion of 
electric dipole and quadrupole transitions in the pre-edge 
can be derived, by measuring XANES spectra for vari- 
ous known orientations (e,k) of the incident beam. 4 In 
addition, the symmetry of the crystallographic sites, that 
host the investigated element, is often a subgroup of the 
cubic group. The number of equivalent sites is given by 
the ratio of the multiplicity of the space group and the 
multiplicity of the point group. The XANES spectrum of 
the cubic crystal is thus the average over the equivalent 
sites of the individual site spectra. Hence, the derivation 
of structural and electronic information for a single site is 
not straightforward, which makes the use of group theory 
and theoretical computations mandatory. 

Among cubic oxides, spinels have attracted consid- 
erable interest for their optical, electronic, mechani- 
cal and magnetic properties. 5-7 In the Earth's interior, 
the formation of silicate spinels has major geophysi- 
cal implications. 8 More specifically, MgAl204 spinels are 
used in a broad range of applications, including opti- 
cally transparent materials, catalyst supports, nuclear 
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FIG. 1: (Color online) Cubic cell of the spinel structure and 
experimental setup. The four equivalent octahedral sites are 
labeled according to the coordinates given in Table I. Site 
distortion has been slightly exagerated for clarity. The sample 
is cut along the (110) plane (red) and rotated along the [110] 
direction, while s and k are fixed in the laboratory system. 
The figure corresponds to the experimental setup taken as 
starting point (a ro t = 0°). For this configuration, the [001] 
axis of the cube is in the vertical plane, perpendicular to 
£ = [010] and k = [100]. 



electronic approach based on density functional theory 
(DFT) and multielectronic methods based on the lig- 
and field multiplct theory (LFM). The monoelectronic 
approach is usually dedicated to the study of delocal- 
ized final states (e.g., the calculation of K-edge spectra) 
but has also provided satisfactory results for the study of 
Ti K pre-edge in Ti0 2 and SrTi0 3 , 16-19 and also for the 
study of Fe K pre-edge in FeS2- 20 The multielectronic ap- 
proach, usually dedicated to the study of localized final 
states (e.g., K pre-edge, L 2j 3 edges of 3d transition ele- 
ments) has been succesfully applied to the case of K pre- 
edge in several systems. 21 ' 22 Our aim is to determine the 
factors (site distortion, electronic interactions) prevail- 
ing at the angular dependence of Cr K pre-edge in spinel 
and to provide a comparison between the monoelectronic 
and multielectronic approaches. We also present a pow- 
erful method, based on symmetry considerations, to re- 
duce the number of calculations needed to reconstruct 
the angular dependence of the cubic crystal from that of 
a single site. The paper is organized as follows. Section 
II is dedicated to the experimental work, including the 
sample description, the X-ray absorption measurements 
and analysis. Section III is devoted to the computational 
work, including the theoretical framework (Sec. Ill A), 
the details of DFT calculations (Sec. Ill B) and of the 
multiplet calculations (Sec. Ill C). Results are presented 
in Sec. IV and discussed in Sec. V. 



waste management and cement castables. 10-12 Cr 3+ of- 
ten substitutes for Al 3+ in MgAl204, which causes a red 
color and makes natural Cr-spinels valuable gemstones. 
Cr 3+ is intentionally added to high-temperature refrac- 
tory materials to improve their thermal and mechani- 
cal properties. 13 In MgAl204 spinel (FdSm space group 
symmetry) , Al 3+ cations occur at octahedral sites, which 
exhibit D%d (or 3 m) symmetry and build chains aligned 
along the six twofold axis of the cubic structure. 14 The 
number of equivalent octahedral sites in the unit cell is 
four, denoted hereafter as sites A, B, C and D, depending 
on their direction of distortion, either [111], [111], [HI] 
or [111], respectively (Fig. 1 and Tab. I). During the Al 
to Cr substitution, the local D$d symmetry is retained 
and the relaxed Cr-site remains centrosymmetric, which 
indicates the absence of Cr 3d-4p mixing. 15 Hence, the 
K pre-edge features arise from pure electric quadrupole 
transitions (Is— *3d) but an experimental evidence of this 
is still lacking. As the Cr-site remains distorted in the 
same direction as for the Al-site, four equivalent relaxed 
sites are available for Cr. Hence, the electric dipole and 
electric quadrupole absorption cross-sections, for a given 
single crystal configuration, are expected to be different 
for a Cr impurity located at site A, B, C or D, since 
their orientations with respect to the incident beam are 
different. 

In this paper, we compare the experimental angular de- 
pendence of the Cr K pre-edge MgAl204:Cr 3+ to those 
obtained by theoretical calculations, combining a mono- 



II. EXPERIMENTS 

A natural gem-quality red spinel single crys- 
tal from Mogok (Burma), with composition 
(Mgo.9 5 Feo.oi)o.96(Al2.o2Cro.oi)2.o304, was investi- 
gated (for details, see Ref. 15). The single crystal was 
cut along the (110) plane (plotted in red on Fig.l) and 
orientated according to the Laue method. 
Cr K-edge (5989 eV) XAS spectra were collected at 
room temperature at beamline ID26 of the European 
Synchrotron Radiation Facility (Grenoble, France). 23 
The energy of the incident radiation was selected using 
a pair of He-cooled Si crystals with (111) orientation. 
The spot size on the sample was approximatively 
250x50 /mi 2 . The orientated sample was placed on a 
rotating holder at 45° with respect to the incident beam, 
and turned around the [110] direction from a rotation 
angle a ro t. The starting configuration (a rot = 0°) 
corresponds to e = [010] and k = [100] (see Fig.l). The 
(&rot — 90°) configuration corresponds to £ = [|,|,-^] 

and k = [-^r^j"^]- F° r this sample cut and this 
experimental setup, the maximum variation effect is 
obtained by substracting the absorption recorded for 
c^rot = 0° from that recorded for a rot — 90°. One 
spectrum was recorded every 15° from a ro t = 0° to 
a rot = 360°, which enables to reconstruct the complete 
angular dependence of the crystal. The absorption was 
measured by a photodiode fluorescence detector. For 
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TABLE I: Coordinates of Cr atom and direction of site distortion for the four equivalent substitutional sites belonging to 
the rhombohedral unit cell. We also give the coordinates of the twelve other sites, obtained from the previous by the three 
translations of the fee lattice (see text and Fig.l). 

site identification direction of site distortion Cr-coordinates in rhombohedral unit cell Cr-coordinates in cubic cell 
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each a ro t angle, ten pre-edge spectra ranging from 5987 
to 5998 eV were recorded with an energy step of 0.05 eV 
and averaged. Two additional scans were recorded 
between 5985 and 6035 eV by step of 0.2 eV, in order to 
merge the pre-edge on the XANES spectrum, and two 
more spectra were recorded between 5950 and 6350 eV 
by step of 0.5 eV, in order to normalize the XANES 
to the K-edge jump far from the edge. Self-absorption 
effects are negligible, because of the low Cr-contcnt of 
the sample. 

III. THEORY 

In this section, we recall the general expressions of the 
electric dipole and quadrupole absorption cross-sections 
for a cubic crystal and for a site with D 3d symmetry 
(Subsec. A). Then, we use the general method described 
in Rcf. 24 to calculate the angular dependence of the cu- 
bic crystal from that of a single site. This framework 
is illustrated in the particular case of spinel. Finally, 
we report the details of the monoelectronic and multi- 
electronic calculations performed for substitutional Cr in 
spinel (Subsec. B and C). 



<7 < ^(4,m)), which transform under rotation like the corre- 
sponding spherical harmonics (Yq, Y™ and Y™). 25 The 
tensor components are functions of Huj, omitted for clar- 
ity in this paper. The symmetries of the crystal restrict 
the possible values of <r D (2,m) and CT^(4,m), as will be 
precised hereafter for the cubic case. 

The electric dipole cross-section for a cubic crystal, 
<J® ub , is isotropic (e.g., it does not depend on the direction 
of the polarization vector) and is equal to &® ub (0,0): 25 

a^(£) = af„ fc (0,0). (2) 

The electric quadrupole cross-section for a cubic crys- 
tal, <rf ub , is expressed, according to group theory (Ap- 
pendix A), as: 

- ^(0,0) + -?L(4^ + 4^ + e 2 fc z 2 -i)al(4,0), 

(3) 

where 0^(0, 0) is the isotropic electric quadrupole cross- 
section, and <jf ub (4, 0) is a purely anisotropic electric 
quadrupole term. The polarization unit vector £ and 
the wave unit vector k have their coordinates expressed 
in the Cartesian reference frame of the cube. 



A. Theoretical framework 

1. Absorption cross-sections for a cubic crystal 

The total absorption cross-section for a crystal (cubic 
or non-cubic), cr, is expressed as: 

<j(£,k) = a D (£)+a Q (£,k) (1) 

where a D is the electric dipole cross-section and a® is 
the electric quadrupole cross-section. The expression 
given above is valid in the absence of coupling between 
the electric dipole and the electric quadrupole terms: 
this condition is fulfilled if the system is either cen- 
trosymmetric or if, at the same time, the system is non- 
magnetic (no net magnetic moment on the absorbing 
ion) and one uses exclusively linear polarization. For 
Cr in MgAl 2 04, the two types of conditions are satis- 
fied. The dipole and quadrupole cross-sections can be 
expressed in function of spherical tensor components, re- 
spectively (ct d (0,0), o- D (2,m)) and (cr Q (0,0), er Q (2,m), 



2. Absorption cross-sections for a site with D 3d symmetry 

For a site with D 3 ^ symmetry, the reference frame is 
chosen consistently with the symmetry operations of the 
point group, i.e. with the z-axis parallel to the C3 axis 
of the D 3 d group. 26 The polarization and the wave unit 

(sin 9 cos <A 
sin6> sincj) 
cos 6 J 

(cos 9 cos 4> cos tp — sin cf) sin tp\ 
cos 9 sin cf> cos ip + cos 4> sin ip \ 
— sin 9 cos ip J 

Hence, 9, which appears in the expression of £ and 
k, is the angle between £ and the C3 axis. 

The electric dipole absorption cross-section in is 
given by: 25 

<M) = <*(0,0) - -±=(3cosH- 1)<(2,0). (4) 
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In order to determine (e) for any experimental 
configuration (e), one needs first to determine cr^ (0,0) 
and cr^ d (2,0), for example by performing calculations 
for at least two independent orientations of i. The 
isotropic term, cr^ (0, 0), can be calculated directly by 
choosing 6 = arccos -j= . 

The electric quadrupole absorption cross-section in 
D 3d is given by: 25 



4- Calculation of the absorption cross-sections for the 
experimental orientations 

The electric dipole isotropic cross-section of the cubic 
crystal, <J^ ub (£), does not depend on the direction of the 
incident polarization vector e. Hence, it will be the same 
for every experimental configuration: 



Vcub( a rot) = V cub 



(0,0). 



(10) 



'D 3d 



(e,k) 



+ — L(35 sin 2 6> cos 2 6 cos 2 -0 



5 sin 2 (9 sin 2 0-4) a"g 3<j (4,0) 



— vlO sin #[(2 cos 6 cos ip ~ 1) cos 6* cos 30 

— (3 cos 2 — 1) sin0 cos ip sin30] cr^ (4,3). (5) 

To determine cr^ 3d (e, k) for any experimental config- 
uration (e,k), one needs first to determine cr^ d (0,0), 

°iL( 2 ' )' a D 3 S^ Q ) and a D 3d ( 4 ' 3 ^ for exam P le b y Per- 
forming calculations for at least four independent orien- 
tations (e,k). 



3. From a single site D 3d to the cubic crystal 

In order to reconstruct the angular dependence of the 
cubic crystal from that of a single site with D^d sym- 
metry, the tensor components have to be averaged over 
the equivalent sites of the cubic cell. For the elec- 
tric dipole cross-section, we need a relation between 
(ag 3d (2,0), <rg 3d (0,0)) and af„ fc (0,0), and for the elec- 
trie quadrupole cross-section, we need a relation be- 
tween (0,0), ^(2,0), ^(4,0), ^(4,3)) and 

<J® ub (4,0)). To do so, we have used the formu- 
las given in Ref. 24, which have been obtained from a 
spherical tensor analysis. This general method uses the 
symmetry operations of the crystal, which exchange the 
equivalent sites of the cubic cell, and is here illustrated 
in the case of spinel. 

The averages over the four equivalent sites are given 
by: 24 



7 D 

cub 



(0,0) 



<,(0,0), 



^(2,0) = 



(6) 
(7) 



For the sample cut and the experimental setup used in 
this study, the expression of the electric quadrupole cross- 
section of the cubic crystal, is given in function of the 
rotation angle a rot by: 



^bKot) = ^(0,0) + ^[-19- 60cos(2a rot ) 

+ 15cos(4« rot )] ( r c Q tlb (4,0). (11) 



The connection between Eqs. 11 and 3 is made following 
the definition of £ and k as functions of a rot (Appendix 
B). Eq. 11 shows that the total angular dependence of 
the cubic crystal is a 7r-periodic function. The fact that 
the rotation axis might not be perfectly aligned with the 
x-ray beam or that the sample might not be perfectly 
homogeneous, could have introduced an additional 2ir 
periodic component. This component would be removed 
from the signal, using a filtering algorithm, based on the 
angular dependence recorded from 0° to 360°. 20 In our 
experiments, this 2tt periodic-component was measured 
to be very small, and no filtering was applied. 

For the present sample cut and experimental setup, 
the maximum variation of the electric quadrupole cross- 
section is expected between a rot = 0° and a rot = 90°. 



• For a r 



0°: 



^ 6 Kot - 0°) = o~q -t=<& 6 (4, 0). (12) 



• For a rot = 90°: 
^Ko* = 90°) - c/q 
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• The isotropic cross-section is : 



<& 6 (4,0). (13) 



1 

15 

= <&6(0,0). 



ol a = ^(8 aZ b (a r ot = 90°) + 7 af ub (a rot = 0°)) 



(14) 



Similarly, we have: 



.24 



7 Q 

cub 



(0,0) 



T D 3d 



(0,0), 



<£*(4,0) = --(7^ m (4,0) + 2V70^ m (4,3))(9) 



(8) 



The dichroic term is : 

vlckro = ^Kot = 90°) - af ub (a rot = 0°) 

15 



2V14 



<&(4,0). 



(15) 
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B. Computational details 

1. Density Functional Theory Calculations 



The computation of the electric dipole and electric 
quadrupole absorption cross-sections were done using 
a first-principles total energy code based on DFT in 
the Local Density Approximation with spin-polarization 
(LSDA). 27 We used periodic boundary conditions, plane 
wave basis set and norm conserving pseudopotentials 28 
in the Kleiman Bylander form. 29 The parameters for the 
pseudopotential generation are given in Ref. 15. 

We started from a host structure of MgAkiCu, which 
is obtained by an ab initio energy minimization calcula- 
tion. In this calculation, the lattice parameter was fixed 
to its experimental value, 30 while the atomic positions 
were allowed to vary to minimize the total energy and 
the interatomic forces. We then relaxed a 2x2x2 rhom- 
boedral supercell containing one Cr atom in substitution 
for Al (i.e., 1 Cr, 31 Al, 16 Mg and 64 O), with the basis 
vectors expressed in a cubic frame. The supercell was 
large enough to avoid interactions between neighboring 
Cr atoms. As the Cr impurity is in its high-spin state, 
the spin multiplet S 2 =| is imposed for the supercell. The 
atomic positions in the supercell were allowed to vary, in 
order to minimize the total energy and the interatomic 
forces. We used a 90 Ry energy cutoff and a single k- 
point sampling in the Brillouin zone. The Cr-site, after 
relaxation, still exhibits a D^,d symmetry, with an inver- 
sion center, one C3 axis and three Gi axis. 

The Cr K-edge absorption cross-section was computed 
using the method described in Refs. 31,32. First, 
we calculated self-consistently the charge density of 
the system, with a Is core- hole on the substitutional 
Cr atom. Then, the all-electron wave functions were 
reconstructed within the projector augmented wave 
framework. 33 The absorption cross-section was com- 
puted as a continued fraction, using a Lanczos basis 
constructed recursively. 34 ' 35 We used a 70 Ry energy 
cutoff for the plane-wave expansion, one k point for the 
self-consistent spin-polarized charge density calculation, 
and a Monkhorst-Pack grid of 3x3x3 fc-points in the 
Brillouin zone for the absorption cross-section calcu- 
lation. For the convolution of the continued fraction, 
we used an energy-dependent broadening parameter 
7, which takes into account the main photoelectron 
damping modes (core-hole lifetime and imaginary part 
of the photoelectron self-energy). The energy-dependent 
7 used in this study is that of Ref. 36. The calculated 
spectrum was then shifted in energy to the experimental 
one: the maximum of absorption is set at 6008.5 eV. The 
absorption edge jump is set to 1, so that experimental 
and calculated spectra for all figures are normalized 
absorption. In such a way, the calculated pre-edge 
could be compared directly to the experimental one. 
As mentioned previously, the four substitutional sites 
will exhibit different spectra for the electric dipole 
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FIG. 2: (Color online) Pre-edge spectra calculated for site 
A, B, C and D in the electric quadrupole approximation. 
The orange lines present the spectra calculated for a rot = 0° 
(e = [010], k = [100]), while the black lines are the spectra 
calculated for a rot = 90° (e = [5, 5,^], k = [-1,-1, -L]). 
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and quadrupole operators, since their orientations are 
different with respect to the incident beam absorption of 
x-rays with given (e,k). The general method to obtain 
the angular dependence measured for the cubic crystal 
is to compute the electric dipole and electric quadrupole 
absorption cross-sections for a Cr impurity lying in each 
of the four trigonally distorted sites A, B, C and D, and 
then to take the average. However, this heavy brute force 
method requires the calculation of four monoelectronic 
potentials with core-hole (after previous associated 
structural relaxation). The number of calculations 
can be drastically reduced if we take advantage of the 
symmetry properties of the crystal, which enables to 
perform the calculations for only one substitutional site 
(site A, with coordinates of (0, j, |) and direction of 
distortion [111]). This method is detailed in Appendix C. 



Figure 2 presents the normalized electric quadrupole 
cross-section calculated for the four equivalent sites, for 
a rot — 0° and a rot — 90°. The spectra calculated for 
ctrot — 0° (orange line) are equal for the four sites. For 
ce ro t — 90°, sites B and C give the same spectra (black 
solid line), as well as sites A and D (black dashed line). 
We observe a slight difference in intensity for the peak at 
5993.2 eV: this is indeed a consequence of the fact that 
sites (A, D) and (B, C) have different orientations with 
respect to the incident beam, and that their symmetry 
differs from Oh- Because the trigonal distortion of the 
octahedra is small in spinel, the anisotropic behaviour 
of the sites is limited for the investigated configurations. 
However, the effect of the trigonal environment can have 
drastic consequences when the distortion is more pro- 
nounced. 
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2. Ligand Field Multiplet Calculations 

In order to extract quantitative information from 
the angular dependence of the pre-edge, we have per- 
formed LFM calculations using the method developed 
by T. Thole in the framework established by Cowan and 
Butler. 26 ' 37 ' 38 In this approach, Cr 3+ is considered as 
an isolated ion embedded in a crystal field potential. 
The band structure of the solid is not taken into ac- 
count, which prevents to calculate transitions to delocal- 
ized (i.e., non-atomic) levels. In other words, the LFM 
approach can be used to calculate K pre-edge spectra, but 
the edge region cannot be computed. Since the Cr-site 
is centrosymmetric, no hybridization is allowed between 
the 3<i-orbitals and the 4p-orbitals of Cr. Hence, the pre- 
edge is described by the transitions from the initial state 
ls 2 3d 3 to the final state ls^d 4 . 

We expose briefly the principles of multiplet calcula- 
tions but details can be found in other references (see 
for example Ref. 39). This approach takes into account 
all the 3d-3d and ls-3rf electronic Coulomb interactions, 
as well as the spin-orbit coupling on every open shell of 
the absorbing atom, and treats its geometrical environ- 
ment through a crystal field potential. In the electric 
quadrupole approximation, the spectrum is calculated as 
the sum of all possible transitions for an electron jumping 
from the Is level toward one 3d level according to: 

<r Q (£,k) = 7T 2 k 2 aftw V— UFle-rk-rinfSiEF-Ej-fvu), 

(16) 

where |I) and |F) are the multielectronic initial and 
final states, of respective energies E/, Ep, and di the 
degeneracy of the initial state. 40 . 

Once the 1 1) and |F) states have been calculated, the 
absolute intensities of the pre-edge spectra are calculated 
in A 2 at T = 300 K. The population of the ground-state 
levels |I) is given by a Boltzmann law. The spectra are 
convoluted by a Lorentzian (with HWHM = 0.54 eV) 
and a Gaussian (with FWHM = 0.85 eV), which respec- 
tively take into account the lifetime of the Is core- hole 
for Cr and the instrumental resolution. Finally, the tran- 
sitions are normalized by the edge jump at the Cr K edge, 
calculated for a Cr atom from Ref. 41 as 4.48 10" 4 A 2 . 
Hence, the calculated spectra can be directly compared 
to the normalized experimental ones. 

The electric quadrupole absorption cross-section was 
calculated for a Cr 3+ ion lying in symmetry, accord- 
ing to the method described above. The crystal-field pa- 
rameters used in the calculation are those derived from 
optical absorption spectroscopy (D q = 0.226 cV, D a = - 
0.036 eV, D T = 0.089 eV). 42 We used the scaling factor of 
the Slater integrals (k= 0.7), related to B and C Racah 
parameters, given in the same reference. The only ad- 
justable parameter is the absolute position in energy. 

As mentioned in Sec. Ill A 2 (Eq. 5), one needs 
first to determine crg 3d (0,0), crg 3d (2,0), erg 3d (4,0) and 



cr§ 3(j (4, 3), in order to determine <J® ub (0, 0) and tr£ b (4, 0) 
using Eqs 8 and 9. This is done by performing four mul- 
tiplet calculations, which provide four independent val- 
ues of the electric quadrupole cross-section (Appendix 
D). Once this first step has been performed, we used 
Eq. 8 and 9 to derive er^ ub (0,0) and 0-^,(4,0). The 
electric quadrupole cross-section of the cubic crystal can 
then be calculated for any experimental configuration 
using Eq. 11. The electric quadrupole cross-section 
for a ro t — 0° and a ro t = 90°, the dichroic and the 
isotropic spectra were determined respectively according 
to Eqs 12-15. 

IV. RESULTS 
A. DFT calculations 

1. Comparison with experiment 



The XANES spectrum, calculated for the cubic 
crystal by first-principles calculations (solid line), is 
shown in Fig. 3 and compared to the experimental 
spectrum (dotted line). As we mentioned above, the 
main absorption edge is due to electric dipole transi- 
tions. Hence, the XANES spectrum does not show any 
angular dependence, except in the pre-edge region. The 
agreement between the experimental and theoretical 
spectra is good, since all the features are reproduced by 
the calculation. A more detailed discussion is reported 
in Ref. 15. The inset of Fig. 3 shows the theoretical 
isotropic XANES spectrum in the pre-edge region (black 
solid line). This spectrum is the sum of the isotropic 

Absorption (arb. units) 




5990 6000 6010 6020 6030 6040 6050 
E(eV) 

FIG. 3: (Color online) Comparison between experimental 
(dotted line) and calculated (solid line) isotropic XANES 
spectra at the Cr K-edge in spinel. The calculation was per- 
formed using DFT-LSDA (see Sec. IIIB1). The inset presents 
the spectra in the pre-edge region. The dashed orange line is 
the calculated electric dipole contribution. 
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FIG. 4: (Color online) Comparison between experimental 
(dotted line) and calculated (solid line) Cr K-pre edge spectra 
in spinel, using DFT-LSDA. The orange lines correspond to 
i =[010], k = [100] (ar O *=0°). The black lines correspond 
to e = [§,§,^], k = [-5,-1,75] ("rot = 90°). The dichroic 
signal (green lines) is the difference between the black and red 
lines, for the experimental (dotted) and the calculated spectra 
(solid), respectively. 



electric dipole (orange dashed line) and the electric 
quadrupole contributions. Our calculations show that 
electric dipole transitions do not contribute to the 
pre-edge, except by a background, which is actually the 
tail of the absorption edge (Is— *p transitions). This is 
a clear confirmation that Cr K pre-edge features are 
due to a pure electric quadrupole contribution. In the 
pre-edge region, the calculated isotropic spectrum is 
in satisfactory agreement with experiment, since the 
two features visible in the pre-edge are reproduced. 
Similar calculations have been successfully performed 
to calculate the K pre-edge for substitutional Cr 3+ 
in corundum and beryl, 36 ' 43 with a good agreement 
between the experimental and theoretical data. This 
shows that a monoelectronic approach can reproduce 
pre-edge features, as can be measured on powder spec- 
tra. However, the position of the theoretical spectrum 
is shifted by about 0.9 eV relative to experiment. 
This shift, which has been already observed in several 
systems, 17 ' 20 ' 36 ' 43 is due the limitation of DFT-LSDA 
in the modelling of electron-hole interaction. In the 
calculation, the effect of the core-hole is to shift the 2>d 
levels to lower energy, with respect to the main edge. 
Unfortunately, this effect is not sufficient to reproduce 
the experimental data, because the core-hole seems to 
be partly screened. 16 This could be improved by taking 
into account the self-energy of the photoelectron. 44 

The experimental and calculated pre-edge spectra for 
the two configurations, which give the maximum dichroic 



signal for the sample cut, are shown in Fig. 4. The num- 
ber of peaks is well reproduced in both cases by the cal- 
culation. For both spectra, the intensity of the first peak 
at about 5990.7 eV is close to 4 % of the absorption edge 
on the experimental data, but underestimated by 25 % 
in the calculation. The relative intensity of the peak at 
about 5992.7 eV is overestimated in the 90° configura- 
tion. Additionally, the energy splitting between the two 
peaks is underestimated by the calculation (1.6 eV vs 
2.0 eV experimentally). The small energy shift of the 
first peak between the two configurations, observed as 
positive in the experimental data, is calculated as nega- 
tive. As a consequence of those several discrepancies, the 
theoretical dichroic signal is not in good agreement with 
the experimental one. Compared to the Ti K pre-edge 
calculations in rutile and SrTi03, using a similar mono- 
electronic approach and reported in several studies, 16-19 
the significant discrepancy observed for Cr in spinel may 
seem at first sight unexpected. However, we underline 
the fact that interlectronic repulsions become crucial for 
localized final states (e.g., for ls~^Sd transitions), and 
that Ti has no d electrons in the systems studied. This 
shows that the electronic interactions on the Cr atom are 
too significant to reproduce quantitatively the angular 
dependence of Cr-spincl in a monoelectronic approach, 
although the average description (i.e., the isotropic spec- 
trum) is satisfactory. 

Nevertheless, the monoelectronic calculation is able to 
reproduce the correct number of peaks. Since this mono- 
electronic approach does not take into account spin-orbit 
coupling and does not fully describe the 3d-3d electronic 
repulsion, a monoelectronic chemical vision of an isolated 
Cr 3+ ion can be applied for the interpretation of the cal- 
culated features. 



2. Assignment of the calculated monoelectronic transitions 
within an atomic picture 

In the following, we shall concentrate on the spectra 
associated to the local symmetry (i.e., calculated for site 
A). In the monoelectronic calculations, the spin multi- 
plet S z =| is imposed for the supercell, since the Cr im- 
purity is in its high-spin state. Cr 3+ has an initial elec- 
tronic configuration (tl g ) 3 (e g )° , which means that Cr 3+ 
is a fully magnetized paramagnetic ion in the calcula- 
tion, while it is paramagnetic in the experiment. Indeed, 
it is not possible to impose the fourfold degenerate S=| 
ground-state in the DFT calculation, that requires non- 
degenerate ground states. In order to assign the transi- 
tions visible in the experimental spectra, we would need 
to calculate the average of the spectra for S z =|, S z —- 
§, S z =^ and S z =-7}. However, it is not possible to do 
the calculation for S z =±i in the Kohn-Sham formal- 
ism, since they are linear combinations of three Slater 
determinants. 45 Nevertheless, the spin-polarized compu- 
tation of the XANES spectrum for S z =| enables to un- 
derstand the origin of the pre-edge features: the contri- 
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FIG. 5: (Color online) Electric quadrupole transitions calcu- 
lated for site A, using spin-polarization. The upper panel (a) 
presents the spectrum calculated for (e = ["^Ji^tjjOL k = [-j- - 
-^=,0]). The middle panel (b) presents the spectrum cal- 
culated for a rot = 0° (e = [010], k = [100]). The lower 
panel (c) presents the spectrum calculated for a ro t = 90° 
(e = [■§■>■§■: ^TS 7 ] ' = ["I "I'T^])- F° r eacn configuration, the 
respective contributions of the two spins are plotted (black 
solid line for spin down, orange solid line for spin up), as well 
as the sum (dashed line). The Fermi level is located approxi- 
mately around 5990 eV. 



bution of the two spins (| and |) can be indeed separated, 
which means that we can deduce whether the 3d-orbitals 
have been reached by a Is electron with spin f or J., and 
this for different expressions of the electric quadrupole 
operator. 

As shown in Fig. 1, although the distortion of the oc- 
tahedra has been slighlty exagerated, the oxygen ligands 
are located approximately along the fourfold axis of the 
cube for all the equivalent sites A, B, C and D. Thus, the 
analysis made for site A provides an assignment, which is 
also valid (mutatis mutandis) for the equivalent sites. For 
a given configuration (e, k), we can easily deduce from 
the expression of the electric quadrupole operator which 
3c?-orbital has been probed in this transition. The inter- 
pretation of the features is possible through group the- 
ory in the monoelectronic approach, using the branching 
rules of Oh 3 D^d (Appendix E). The d orbitals belong 
to the t2 g {Oh) and e g {Oh) irreducible representations 
within octahedral symmetry. When lowering the sym- 
metry to Z?3d, the t2g{Oh) irreducible representation is 
split into the eg{D^d) and aig^D^d) irreducible represen- 
tations. To indicate that they come from t2 g (Oh), they 
will be written as e g (t2 g ) and cti ff (i2g)- The e g (Oh) irre- 
ducible representation becomes the e g {D^d) irreducible 
representation, designed hereafter as e g (e g ). 

The normalized electric quadrupole cross-sections cal- 
culated for site A are shown in Fig. 5 for three different 
configurations (e, k). For a better understanding of 
the structures, the electric quadrupole transitions to 
both occupied and empty states are represented. For 

& = [TI'Tl' ]' fc = [7f'-7f'°]) ( Fi s- 5a )> the electric 

quadrupole operator is expressed as O a = -^(x 2 — y 2 ), 
which enables to probe the 3d electronic density in the 
x 2 -]/ 2 direction, i.e., along Cr-0 bonds: the orbitals 
probed are the e s (e g ), which are empty for spin f and 
|, since they are coming from the e g (Oh) levels. Fig. 5a 
shows that, indeed, two peaks are obtained at 5991.6 eV 
and 5993.2 eV, above the Fermi level. Below the Fermi 
level at 5990 eV, a broad structure is observed between 
5982 eV and 5990 eV, which corresponds to e g states 
hybridized with the p-orbitals of the oxygens. 
For a rot = 0° (e = [010], k=[100]), one single peak is 
obtained in the empty states at 5991.6 eV for spin J, (see 
Fig. 5b, black line). For this orientation, the electric 
quadrupole operator, expressed as Ob = xy, enables to 
probe the d electronic density, projected on Cr, in the 
xy direction, i.e., between the Cr-0 bonds. The e g (t2 g ) 
and ciig(t2g) orbitals having a component along xy, as 
indicated by their expressions in Appendix D (Eq. 44), 
they are probed in the transition. As these states coming 
from the t\ g (Oh) are fully occupied, they can be reached 
only by a photoelectron with spin |. This is indeed 
consistent with our results. The splitting between 

and a\ g is not visible, which is an indication of a small 
trigonal distortion for the Cr-site in spinel. Below the 
Fermi level, a broad structure with an intense peak at 
5988.7 eV is observed. The intense peak corresponds to 
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the occupied e g (t2 g ) levels. To interpret the origin of 
the broad structure, we have to remind that the e g (t2 g ) 
states can hybridize with the e g (e g ) levels, since they 
belong to the same irreducible representation in D^d- 
As mentioned previously, the hybridization of the mixed 
e g states with the p-orbitals of the oxygens gives rise to 
the structures visible below 5988 eV. 
For a rot = 90° (£ = [§,§,^], k = [-|,-|,^]) (Fig. 5c), 
two peaks are obtained above the Fermi level. For this 
orientation, the electric quadrupole operator is expressed 
as O c — 3z ^ r — which enables to probe the 3d elec- 
tronic density both in the xy and 2>z 2 — r 2 directions. 
For the 3z 2 — r 2 component, the levels probed are the 
e g (e g ), as for the first orientation studied (Fig. 5a). For 
the xy component, the levels probed are the e 9 (<2g) and 
Q>ig(t2 g ), as for the second orientation (Fig. 5b). Fig. 5c 
shows that the spectrum is a close combination of the 
transitions visible on the two previous spectra (Fig. 5a 
and b), and the assignment of the structures is made 
clear from the two previous cases. The position of the 
4g( e s) peak is close to that of the e g (e g ) at 5991.6 eV. 
The energy difference (1.6 eV) between 4 ff ( e s) an d 
e g (e g ) gives an idea of the t2 g (Oh)-e g (Oh) splitting 
due to the crystal field. This can be compared to the 
experimental crystal-field splitting (2.26 eV), derived 
from optical absorption spectroscopy, but one should 
keep in mind that the crystal field splitting in the 
monoelectronic picture is associated with spin | levels. 
For the configuration O c , the 3z 2 —r 2 component enables 
to probe the e g (Oh) states, as x 2 —y 2 . Considering the 
normalization factors in the expression of the d orbitals, 
the magnitude of the transition operator along 3z 2 — r 2 
is y/3 times bigger than the magnitude of the transition 
operator along x 2 — y 2 , which is 2 times bigger than the 
magnitude of the transition operator along xy. O c thus 
appears as a linear combination of the two operators 
O a and Ob, with respective weights of and |. If no 
coupling occurs between O a and Ob when calculating 
the square matrix elements |(/|e • rk • r|z)| 2 , the third 
cross-section (c) should be the linear combination of 
the two cross-sections (a) and (b) obtained for O a and 
Ob, with respective weights of | and j. However, 
the linear combination of the two cross-sections (not 
shown in Fig. 5) and the cross-section obtained for 
the linear combination of the transition operators are 
slightly different, which indicates a small interference 



between the xy and Zz A 



(or x — y ) compo- 



nents. The interference is a clear evidence of the e g {e g ) 
and e^(i2 9 ) hybridization due to the D^d local symmetry. 



B. LFM calculations 

1. Comparison with experiment 

For the two experimental configurations, Fig. 6 
presents the experimental Cr K pre-edge spectra (dot- 
ted line), the theoretical spectra obtained by LFM cal- 
culations (solid line) and the corresponding dichroic sig- 
nals. The calculated pre-edges have been obtained for 
the cubic crystal from a calculation performed for a sin- 
gle site with D$d symmetry. For each configuration, the 
shape of the spectrum is well reproduced by the calcula- 
tion. In the experimental data, the position of the first 
peak is shifted by approximately +0.15 eV for a rot = 0°, 
compared to that in the a ro t = 90° configuration. This 
relative shift is also well reproduced in the calculated 
spectra. For a ro t = 90°, the relative intensity of the two 
peaks is in good agreement with the experimental data. 
The shape of the dichroic signal is well reproduced by 
the calculation: in fact, the x-ray linear dichroism of the 
crystal is well described in the multiplet approach, sug- 
gesting that the calculation includes the necessary mul- 
tielectronic interactions on the Cr atom. We recall that 
the crystal-field parameters used in the calculation are 
those obtained from optical absorption spectroscopy (see 
Appendix F & G for the correspondance between the ex- 
perimental crystal-field parameters and the parameters 
used in the multiplet calculations). However, the inten- 
sity of the dichroic signal is overestimated by 20 % in 
the calculation. The first reason for this overestimation 
is that the calculated spectra have been normalized by 
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FIG. 6: (Color online) Comparison between experimental 
(dotted line) and calculated (solid line) Cr K pre-edge spectra 
in spinel, using Ligand Field Multiplet theory in D^d symme- 
try. The orange lines correspond to £ — [010], k = [100] 
(ctrot = 0°). The black lines correspond to £ — [^,^,^], 
k = [-§,-§1775] {etrot = 90°). The dichroic signal is the differ- 
ence between the black and orange lines. 
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the edge jump at the Cr K edge, which was calculated 
for an isolated Cr without considering the influence of 
the crystal structure according to Ref. 41. This can ac- 
count for a few percent in the discrepancy. Another few 
percent possibly lie in the normalization of the experi- 
mental data, since we used the average of two spectra, 
which were recorded between 5950 eV and 6350 eV with 
a rather large energy step (0.5 eV). This can introduce 
limited noise and thus uncertainty on the normalization. 
A third source of error is that the crystal-field parame- 
ters used in the calculation might be slightly different in 
the excited state than in the ground-state, because of the 
influence of the core-hole: for example, if D q is increased 
by 2 % in the excited state, the intensity of the first peak 
in the dichroic signal decreases by 14 %. The shape and 
intensity of the calculated dichroism are quite sensitive 
to the crystal-field parameters used in the excited state. 

Nevertheless, despite this slight intensity mismatch 
with the experimental data, the angular dependence of 
the crystal is well reproduced by the calculation, which 
means that the multielectronic approach takes into ac- 
count the necessary interactions. Since isotropic and 
dichroic calculated spectra fit well with experiment, the 
analysis of the calculation is very likely to yield valuable 
insight into the origin of the experimental transitions in 
the pre-edge region. In the following, we shall investigate 
the influence of the different terms in the Hamiltonian 
taken into account in the LFM approach (trigonal dis- 
tortion, fourfold degeneracy of the ground state trigonal 
S=| ( 4 A2g), spin-orbit coupling on the 3d levels, 3d-3d 
or 1s-3g! Coulomb repulsion) on the angular dependence. 



2. Influence of trigonal distortion on dichroism 

In this paragraph, we investigate the influence of 
the trigonal distortion on the angular dependence. The 
isotropic and dichroic spectra in Oh symmetry have been 
obtained by setting the trigonal distortion of the crys- 
tal field to zero. They are compared to those calculated 
in D^d symmetry using the distortion parameters given 
by optical absorption spectroscopy (D a = - 0.036 eV, 
D T = 0.089 eV). As shown in Fig. 7, the difference be- 
tween the calculations performed in Oh and D^d sym- 
metries (orange solid line and black dotted line, respec- 
tively) is weak, since the isotropic and dichroic signals 
have similar shape and intensity. This result is consis- 
tent with the small values of the parameters D a and D T , 
which quantify the trigonal distortion of the Cr-site in 
spinel. This means that, provided that trigonal distor- 
tion is limited, the calculation of pre-edge spectra could 
have been performed for a single site with Oh symmetry 
(see Appendix A2 for simplified formula). This is also in 
line with the monoelectronic calculation, for which the 
splitting between e|(t2 9 ) and a| g (t 2ff ) could not be re- 
solved in the calculated spectra. 

We have investigated the effect of the intensity of the 
trigonal distortion on the calculated spectra by choosing 
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FIG. 7: (Color online) Comparison between calculated Cr 
K pre-edge spectra in spinel in D^d symmetry for increasing 
trigonal distortion. The orange solid line, labeled Oh corre- 
sponds to a perfect Oh symmetry {D a = 0, D T = 0). The 
black dotted line,labeled D^d corresponds to {Da = -0.036 
eV, D T = 0.089 eV), which are the distortion parameters 
obtained from Optical Absorption Spectroscopy. 42 The black 
solid line, labeled D$d-2, was obtained for a doubled distor- 
tion parameters {D a = -0.072 eV, D T = 0.178 eV), and the 
orange dotted line, labeled D^d-3 for tripled distortion pa- 
rameters {D a = -0.108 eV, D T = 0.267 eV) 



two other sets of the distortion parameters. In Fig. 7, 
the spectra labeled D^-i {i = 2,3) are calculated with 
the set of distortion parameters (i x D a , i x D T ), for 
{D a = - 0.036 eV, D T = 0.089 eV). The crystal-field pa- 
rameters used are the same in the ground- and excited- 
state. As seen in Fig. 7, the intensity of the isotropic 
spectra {D$d, D^d-2 and D^d-S) are almost identical, in- 
dicating that the isotropic signal is not sensitive to site 
distortion. The intensity of the maximum at 5990.75 eV 
remains close to 2.5 % of the electric dipole edge jump. 
On the contrary, the shape and intensity of the linear 
dichroic signal is highly sensitive to the trigonal distor- 
tion. The intensity of the first feature at 5990.25 eV is 
lowered when the distortion is increased: one observes a 
20 % decrease when the distortion parameters are dou- 
bled (signal labeled D 3 d~2 in Fig. 7), and a 50 % decrease 
when the distortion parameters are tripled (signal labeled 
Dzd-3). It should be noticed that, for our parameter sets, 
the increase of site distortion is accompanied by a rather 
counter-intuitive decrease of the intensity of the linear 
dichroic signal, thus indicating the relevance of the the- 
oretical developments performed within this paper. This 
means that site distortion has to be carefully taken into 
account when calculations are performed to mimic the 
angular dependence of the pre-edge. In that case, the 
calculation for a single site with D%d symmetry should 
follow the method described in Sec. III. 
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3. Influence of ground-state degeneracy, spin-orbit coupling 
and interelectronic repulsion on dichroism 

Beyond the site symmetry distortion, the other ingre- 
dients of the calculation are the fourfold degeneracy of 
the S=| ground state, the 3d-3d Coulomb repulsion, the 
ls-Sd Coulomb repulsion and the 3d spin-orbit coupling. 
We shall check the influence of these different parame- 
ters. We have performed multiplet calculations restrict- 
ing the ground-state to the non-degenerate S z =| state of 
the S=| multiplet. The calculated electric quadrupole 
transitions are almost identical to those with the four- 
fold ground state. Differences are below 0.1% of the 
maximum intensity of the isotropic electric quadrupole 
spectrum. This clearly indicates that the procedure fol- 
lowed in monolectronic calculations to take into account 
the spin degeneracy is sound and appropriate. The ra- 
dial integrals for Coulomb interaction and spin-orbit cou- 
pling are calculated by relativistic Hartree-Fock atomic 
calculations. One finds: the ls-3d exchange Slater in- 
tegral G\ s 3d = 0.052 cV, the direct 3d-3d Slater inte- 
grals Fg d 3d = 10.78 eV and F^ d 3d = 6.75 eV, and the 
3d spin-orbit coupling ( 3 d — 0.035 eV. G\ s 3d is small 
compared to F^ d 3d and F* d 3d . Using G\ s3d = in the 
multiplet calculation, we found almost no difference with 
the isotropic and dichroic signal calculated with the ab 
initio atomic value of G\ s 3d . By calculating the dichroic 
signal with ( 3d ~ 0, we found a small difference concern- 
ing the intensity of isotropic and dichroic signals, when 
compared to the associated spectra with = 0.035 eV. 
The maximum relative difference is less than a few per- 
cent (2 %) of the feature intensity. The observed small 
dependence of the pre-edge features with G\ 8 3d and C,3d 
is in line with results obtained at the Fe K pre-edge. 21 
We also performed calculations setting F 3d 3d and F 3d 3d 
to zero, and we observed that the isotropic and dichroic 
calculated spectra (not shown) were in complete disagree- 
ment with experimental data. This clearly indicates that 
the direct Slater integrals on the 3d shell, and thus the 
multiclectronic 3d-3d Coulomb interactions, are the es- 
sential ingredients governing the shape of the isotropic as 
well as the dichroic signals. 

From the preceding analysis, we can unambiguously 
determine the parameters governing the shape and 
intensities of the pre-edge features. Spin-orbit coupling 
on the 3d orbitals, ground-state degeneracy and ls-3<i 
Coulomb repulsions have only limited impact on the 
calculated LFM isotropic spectra. This explains the 
reasonable agreement between calculation and exper- 
iment for isotropic pre-edge in the DFT formulation, 
where the two first previous ingredients arc missing, 
and where the ls-3d Coulomb repulsion is taken into 
account in an approximate way. The D^d distortion has 
almost no influence on the isotropic pre-edge but can 
have a large one on the dichroic signal. In the case of 
Cr in spinel, the trigonal distortion is such small that 
it does not provide detectable features on the dichroic 



signal. The major ingredient for the interpretation of 
the Cr pre-edge features is 3rf-3<i Coulomb repulsion. 
This effect is highly multielectronic and complicates the 
simple interpretation provided by the monoelectronic 
scheme. This ingredient is mandatory to get correct 
intensities and energies for both isotropic and dichroic 
signals. 



V. DISCUSSION AND CONCLUSION 

On the one hand, monoelectronic calculations allow 
to make contact between electric dipole and electric 
quadrupole calculations. They show that electric dipole 
transitions do not contribute to the features visible in the 
pre-edge and they provide a clear vision of the assign- 
ment of the \s-3d transitions occuring in the pre-edge. 
However, they are unable to reproduce quantitatively the 
linear dichroism in cubic crystals, since the interelec- 
tronic repulsion on the 3d levels of the Cr ion cannot 
fully be described in the LSDA framework. On the other 
hand, multielectronic calculations well reproduce the an- 
gular dependence of the pre-edge in cubic crystals, as well 
as the isotropic spectrum, with no adjusted parameters. 
However, in this approach, the main absorption edge, 
associated to electric dipole transitions, cannot be repro- 
duced since the band structure (or at least the electronic 
structure of a large enough cluster around the absorbing 
atom) is not taken into account. The agreement between 
experiment and multiplet calculations indicates that the 
assignment of the transitions is no more straightforward, 
as could have been expected from a more simple atomic 
monoelectronic picture. Hence, the two approaches are 
highly complementary. 

From this monoelcctronic-multiclectronic combined 
approach, our first finding is that the 3d-3d electronic 
repulsions and the crystal field are the main interactions 
prevailing at the K pre-edge of Cr in spinel. The mul- 
tiplet approach seems mandatory to describe quantita- 
tively the K pre-edge of 3d transition ions, and more 
generally the K-edge spectra of elements for which elec- 
tronic correlations are significant. The effect of the 3d 
spin-orbit coupling and of the ls-3<i Coulomb repulsion 
are very weak: this explains that the monoelectronic ap- 
proach (which does not fully take into account these in- 
teractions) can provide a satisfactory simulation of the 
isotropic spectrum. 

Our second finding concerns the maximum proportion 
of electric quadrupole transitions in the Cr K pre-edge 
that can be estimated, with respect to the edge jump 
(here, normalized to 1): the intensity of the largest peak 
(5990.75 eV) on the pre-edge isotropic spectrum is less 
than 2.5 % of the edge jump. From the monoelectronic 
calculations, we estimate that the non-structured slope 
from electric dipole origin contributes to about 0.9 % of 
the edge jump at 5990.75 eV. Thus, the total intensity 
of the largest pre-edge feature does not exceed 3.5 %. 
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We can conclude that, if the pre-edge features are more 
intense than 4 % of the edge jump, pure quadrupolc 
transitions alone cannot explain the origin of the struc- 
tures. It gives a strong limitation to the often encoun- 
tered idea that electric quadrupole transitions could ex- 
plain large pre-edge features. This result, which is consis- 
tent with previous studies on Fe 2+ and Fe 3+ in minerals 
and glasses, 21 ' 22,46 can probably be extended to the other 
3d transition ions. 

Our final finding concerns the relation between the 
spectral feature of the pre-edge with the local site 
distortion of the absorbing ion. The effect of the trigonal 
distortion does not affect significantly the pre-edge 
isotropic spectrum. This is a general trend already 
observed for other related spectra (electric dipole tran- 
sitions for K-edge and L 2 ,3 edges of 3d elements). On 
the contrary, the dichroic signal is much more sensitive. 
This indicates the possibility to obtain quantitative 
information on site distortion from the linear dichroic 
dependence of the pre-edge feature. This can only be 



recorded if angular dependent measurements on single 
crystals are performed to yield the full dependence 
of the absorption signal. The connection between 
site distortion and linear dichroism is then made by 
simulations in the LFM method within the geometrical 
analysis developed throughout this paper. 
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VI. APPENDIX 

A. Proof of Equation 3 

1. General expression of the electric quadrupole 
cross-section for a cubic crystal 

We start from the defining formula of the electric 
quadrupole cross-section for linearly polarized x-rays: 

cr Q (f,k) = Tr 2 k 2 ahuj^2\(f\(s-rk-r\g)\ 2 S(E g + hLu-E f ), 
f 

(17) 

ijlra 

where : 

Vijlm = T: 2 k 2 ahujY,}{9V l n\f){.f\rjr m \g)5{E l + hhj-E f ) 1 

£ and k arc the polarization and wave unit vectors, 
respectively. 

To calculate the form of this sum when the sample has 
a symmetry group G, we use the fact that the absorp- 
tion cross-section is invariant by any symmetry opera- 
tion that acts on both the sample variables (Jijim and 
the x-ray variables £, k. Therefore, the cross-section 
is left invariant by the crystal symmetries applied to 
the x-ray variables. In other words, for any operation 
S of the symmetry group G of the sample, we have 
a(£, k) = a(S(£), 5(k)). Therefore, if G is the symmetry 
group of the sample, or a subgroup of it, we can write 



^(£,k) 



i 

IGI 



5>(S(S),S(k)), 

s 



(18) 



where \G\ denotes the number of elements of G. 



We rewrite Eq. 17 as 

+ Si^j £ i £ jkjkiCFijji + e i k t ana + R 



^£ jk^kjO \; 



(19) 



Eq. 19 defines term R. The term R is the sum of the 
terms that are not (i = j and I = m) or (I = i and 
to = j) or (ttt, = i and I = j). So R is a sum of 4 
terms of the type (i 7^ j and i 7^ I and i 7^ 777) plus 
the three cyclic permutations of (i,j,l,m) and 4 terms 
of the type (i 7^ j and j 7^ I and I 7^ i) plus the three 
cyclic permutations of (i, j, l,m). We want to prove Eq. 3 
that is valid for a cubic crystal in a reference frame such 
that x, y, z axis are taken along the fourfold symmetry 
of the cubic crystal. We first show that if the sample has 
three perpendicular mirror planes, the term R is zero. 
Consider the term (i 7^ j and i 7^ I and i 7^ 777) and take 
the symmetry (S(si) — — £j) the other variables j, k, I are 
different from i, so the symmetry leaves them invariant 
and changes only one sign. Therefore, using Eq. 18, this 
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term disappears from R. The same is true for the three 
cyclic permutations. Consider now the term (i ^ j and 
j I and I ^ i). Since the values of the indices is 1, 2 or 
3, one of the three indices i,j,k is different from the other 
two and from m. So one index is again different from the 
other ones and the same reasoning can be applied to show 
that the corresponding term vanishes. This holds also for 
the three cyclic permutations and we have shown that, 
when there are three perpendicular symmetry planes, the 
absorption cross-section is 



+ (20) 



The group Oh has a subgroup made by the six permu- 
tations of (x,y, z). An average over this subgroup gives 
the following result 



cub 



{£, k) = £^,- elk) A + £. . . e l e j k i k ] B 



E^kfC, 



(21) 



where 



A 
B 

C 



3 ' 



To complete the proof, we use the fact that 



^ ^ SiSjkikj 



(f.k) 2 -^ e ^, 



and the identities e-£ = l,k-k=l and £ ■ k = to get 
<7l(£,k) = yL + ]T £ 2 fc 2 (C-.4-i?). (22) 

To compare this result with the expansion over spher- 
ical tensors, we need to determine the isotropic contri- 
bution Uq , which is obtained as the average of af ub (£, k) 
over angles. We write a(e, k) in terms of 9, 0, ip as in 
Sec. Ill A 2. 

Thus, 



a^ b (f,k) = A+ c -^- i? sin 2 6>(7cos 2 6'cos 2 V 
+ cos 40 cos 2 cos 2 ip) + 2 sin 2 20 sin 2 ^ — cos sin 40 sin 20) , 

and the average over all directions is 

*0 = 8^ Jo" Sin ^ Jo ' ^ Jo ' d ^£ 6 ( f ' fc ) 



Therefore, 



and 



k) = + ( £ 2 fc 2 + £ 2 fc 2 + £ 2 fc 2 - (23) 

where erf = (C — A — B) and k x , k y , k z are the 
coordinates of k. 

This expression is valid for any cubic crystal, providing 
that the reference frame is such that the x, y, z axis 
are taken along the fourfold symmetry axis of the cubic 
crystal. 



2. Absorption cross-sections in Oh symmetry 

For a site with Oh symmetry, the orthonormal refer- 
ence frame chosen is that of the cubic crystal. The z-axis 
of the reference frame is parallel to the fourfold axis of 
the cube. The angle 9 is thus the angle between the po- 
larization vector and the z-axis of the cube. 

The cross-section calculated for a single site in Oh 
is equal to the cross-section of the cubic crystal, since 
a perfect octahedron and the cube have the same 
symmetry operations: 

ag fc (e,k) = ^(e.k) and ag h (£) = a° b (£). 

For a single site with Oh symmetry, the expression of 
the electric dipole cross-section is very simple: 



ag h (£) = a D (0,0), 



(24) 



where <r D (0,0) is the isotropic electric dipole cross- 
section. 



For the cubic crystal, one obtains: 



(25) 



The electric quadrupole absorption cross-section for a 
site with Oh symmetry is given by: 



crS (e, k) = ctS (0, 0) + [35 sin 2 9 cos 2 9 cos 2 ip 



1 



= A + 



C-A-B 



14 

+ 5 sin 2 9 sin 2 0-4 

+ 5 sin 2 9 (cos 2 9 cos 2 ?/; cos4 — sin 2 cos 40) 
-2 cos 9 simp cos0 sin 40] <rg h (4,0), (26) 

where <J® uh (0, 0) is the isotropic electric quadrupole 
cross-section, and c^ b (4,0) a purely anisotropic electric 
quadrupole term. 
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Using ag h (e,k) = <7« b (e,k) and (9 = 
?), Eq. 26 is equivalent to Eq. 23 with: 



^ =(7^(0,0) =(7^(0,0), 



Q 



(27) 



n^ (4 ' 0)= t^^ (4 ' 0) - (28) 



Eq. 23 can be rewritten as: 



<&6(e,k) = <£*(0,0) 

+ -^( £ ^+e^+e^-i)al(4,0). 

(29) 

This is the proof of Eq. 3. 



— rfi 



B. Expression of £(a ro t) and k (a ro t) 

For the sample cut and the experimental setup used in 
this study, we have: 



i{otrot) = 



k(a rot ) 



' 1 — COS OL r ot ' 

2 

1+cos a r oi 
2 

Sill OLrot 

V2 

-1 — COS Q. ro 

2 

- 1 + COS OL ro 

2 

sin a ro t 

v5 



C. Symmetry adapted method used in 
monoelectronic calculations 

Firstly, the absorption cross-section was calculated in 
the electric dipole approximation, in order to derive 
o~® ub (0, 0), the isotropic electric dipole cross-section. This 
term can be obtained from a single calculation of the elec- 
tric dipole absorption cross-section performed for site A. 
The expression of the electric dipole cross-section in D 3 d, 
a®, is given by Eq. 4, where 9 is the angle between the 
polarization vector and the C3 axis. The program we 
used for the ab initio calculations calculates the average 
of a® (£) for £ along the x-, y- and z-axis of the cubic 
frame. The angle 9 between the [111] direction (parallel 
to the C3 axis of the site) and each of these three direc- 
tions 

vector £ taken along the x 



is arccos-^. This implies that, for a polarization 



y- or 2-axis, 



°?uM = = < d (ci, 0) = Oo, 0). (30) 

This means that the average value calculated by the pro- 
gram is directly equal to <J® ub (0, 0). Hence, cr^ ub (0, 0) was 
obtained from a single calculation performed at site A. 



Secondly, the calculation was performed in the electric 
quadrupole approximation, in order to derive <j® ub (0,0) 



and o-® ub {&, 0). Once these two terms are determined, 
we will be able to calculate the electric quadrupole ab- 
sorption cross-section for the cubic crystal, according to 
Eq. 3, for any (e,k) configuration. We used a symmetry 
adapted method to determine AO, 0) and cr® b (4, 0) in 
order to reduce the number of calculations: this way, it 
is possible to consider one single substitutional site (site 
A) and take advantage of the symmetry properties of the 
crystal. This method allows to save significant computa- 
tional time, since we perform only two self-consistent cal- 
culations (instead of eight with the brute force method): 
one calculation to do the structural relaxation of the sys- 
tem (substituted at site A) , and a second one to calculate 
the charge density with a core-hole on Cr. As mentioned 
in Ref. 24, assuming that we have calculated the spec- 
trum for a given site X, it is possible to obtain the spec- 
trum for any site Y equivalent to X by calculating the 
spectrum of site X for a rotated x-ray beam. More pre- 
cisely, in the case of electric quadrupole transitions, if site 
Y is the image of site X by a rotation R, the spectrum of 
site Y for a configuration (f,k) is equal to the spectrum 
of site X for the rotated configuration (R _1 (e),R _1 (k)). 
As in Ref. 24, let us consider the site with reduced co- 
ordinates (0, 1/4, 3/4), with is a representative of site 
A. By applying the three rotations about the z-axis of 
the crystal through angles of n/2, w, 37r/2, we obtain 
the positions of three sites, which are respectively rep- 
resentative of sites B, D and C. The three rotations will 
denoted as R^, Rx and R 3t / 2 - This implies that 



*g(S,k) 
ag(f, k) 



a Q A {R-] 2 {£),R-}^)), 

ai(R^ /2 (£),R-^)), 

ai(R-\£),R-\k)). 



(31) 
(32) 
(33) 



We shall now apply the previous equations to the case 
of the two experimental configurations, a rot — 0° and 
a ro t — 90°, in order to see if additional simplifications 
can be found. To do so, we need to consider the ex- 
pression of the electric quadrupole transition operator Q: 



Q = -£ ■ rk • r. 



(34) 



Note that, in text, the assignment of the calculated^mo- 
noelectronic transitions is discussed using O = £ ■ rk • r. 
We have to consider the absolute value of O (or Q) 
because the cross-section does not depend on the sign of 

6. 

For a rot = 0° (e = [010], k = [100]): 



\£ ■ rk • r| = |i?; / 1 2 (e) • rR-) 2 (k) ■ r| = \R~\£) ■ rR-^k) ■ r| 
= |i?- 7r 1 /2 (£)-ri?; / 1 2 (k)-r|= a ;y. 
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TABLE II: Expression of the electric quadrupole cross-section calculated in D^d group for the four independent orientations 
used to derive erg d (0, 0), erg d (2, 0), crg 3d (4, 0) and trg 3d (4, 3) (see text). The coordinates of £ and k are given in the reference 
frame chosen for Dzd, with the 2-direction parallel to the C3 axis. 

label 8 rp £ k qg M (e,k) 

al arccos(^) f (-^, ^(0,0)4-43^(4,0)4-^1^(4,3) 

«« —(TS + Te) (^A^) (^A^) < d (0,0)- V ^<j2,0)- i i^< d (4,0)-^g 3d (4,3) 

s5 | f f (0,1,0) (-1,0,0) ag 3d (0,0) + V / fag 3d (2,0)+ ^ag M (4,0) 

^ ¥ f - (O'Tf-75) (°'7f'7f) < d m-^ 4 a%j2,0)+ ^< d (4,0) 



Hence, we have: 

^(e,k) = ^( J R; / 1 2 (e)^; / 1 2 (k)) 

= ^( J R; 1 (f), J R; 1 (k)) 



= ^(^ 1 /2(^)^ 37r 1 / 2 (k))- 



This means that: 



a Q A (£, k) - (£, k) - (£, k) - 0% (£, k) . (35) 

For a r0 ( = 0°, we thus need to perform one single calcu- 
lation of the electric quadrupole cross-section, since 



For a rot = 90° (f - [±,±,^] and k = [-\-\,^]): 



(36) 



2' 



\£ rk r| = lii-^e) • rij^k) • r| = |z 2 /2 - (a; - j/) 2 /4| 
and 

l^'/V*) ' ^(k) • r| = \R£ /2 {§) ■ ri?;; 2 (k) • r| 
= \z 2 /2-{x + yf/A\. 

This means that: 

a«(£,k) = <rg(S,k), (37) 
<7f(e,k)=ag(£,k). (38) 

For a rot = 90° , we thus need to perform two calculations 
of the electric quadrupole cross-sections since 



cub 



(£,k) 



<r«(£,k)+<rg(£,k) 



(39) 



As mentioned previously, instead of doing the calcula- 
tion for the two sites A and C, it is more convenient to 
compute the spectrum for site A, for the two orientations 
(e,k) and (R~L (e),R~L (k)). This corresponds to (e 

= tl'I'T/j]' k = ["^'"I'T/f]) and (® = ["I'l'Tj]' k = 
[I'-i'Tf])' respectively. 



D. Method used to perform the multiplet 
calculations 

As mentioned in Sec. Ill A 2 (Eq. 5), one needs 
first to determine (7^(0,0), (7^(2,0), crg 3d (4,0) and 

a% 3d (4, 3), in order to determine (7^,(0, 0) and af ub (4, 0) 
using Eqs 8 and 9. This is done by performing four multi- 
plet calculations, which provide four independent values 
of the electric^quadrupole cross-section, si, s2, s3 and s4, 
where e and k are defined in Table II. The components 

^ M (°>°)> ct L( 2 '°)' ct L( 4 >°) and ff L( 4 > 3 ) are ob - 

tained by a combination of si, s2, s3 and s4, according 
to: 



^(4,0) 



12 4 2 

- sl + - s2+— s3+— si 

5 5 15 15 



2Vu Ay/u 2VTI mVu 



35 



35 



2.1^2 

— si s2 — 

\/5 2V5 3^ 



105 

s3- 



105 



(40) 



These equations have been obtained by inversing the sys- 
tem of equations, which give the expressions of si, s2, s3 
and s4, in function of 0^(0,0), ^(2,0), ^(4,0) 
and crg 3d (4,3) (Table II). 

Once this first step has been performed, Eq. 8 and 9 
are used to derive <jf ub (0, 0) and crf ub (A, 0). The electric 
quadrupole cross-section of the cubic crystal can then 
be calculated for any experimental configuration using 
Eq. 11. 



E. Expression of the d-eigenstates in Dzd 

The d-eigenstates in point group are determined 
by the branching rules of the irreducible representation 
2+(0 3 ) in D 3d . In order to get the complete eigenstates, 
we must consider the 3 D D D 3c i D C 3i subduc- 
tion. To simplify the notation, we will make no use of 
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TABLE III: Branching rules for k 
D 3 D C3 subduction 



2 using the SO3 D O D 



SOS 


-> O -> 


£>3 




c 3 


\k(S0 3 )p(0)a(D- A )\(C 3 )) 


2 




1 




1 


M*0> 


2 


-c 1 -> 


1 




-1 


|c-(*a)> 


2 












|ai(*a)> 


2 


-> 2 -> 


1 




1 


|e+(e)> 


2 


-f 2 -> 


1 




-1 


|e-(e)) 



I JM) partners (Rcf. 26 p. 527), we obtain: 
|2200) 4 = -|20) 

|2222) 4 = --^|22) - -J=|2-2) 



|21H>4 = 
|211-1)4 



|21) 
|2-1) 



|2122) 4 = -L\22) - -^|2-2) 



(42) 



the parity (± or g/u) in the rest of the appendix. Hence, 
we will use the subduction SO3 D O D D 3 D C3. In 
the following, the irreducible representations are labeled 
according to Ref. 26. For example, in O group, the irre- 
ducible representation 1 designates T 2 in Schonflies no- 
tation, while 2 designates E. The complete eigenstates 
are written as \k(S03)p(0)a(D3)X(C3)}3, where A is the 
irreducible representation of C3 subgroup, coming from 
the k irreducible representation of SO3, that becomes p 
in 0, a in D3 and A in C3 . The branching rules for k— 2 
are given in Tab. III. 

We recall that in D 3 (or C3), the reference frame 
chosen is not the one used in O. Therefore, we 
need to express the \k(S0 3 )p(0)a(D 3 )X(C 3 )) 3 
in function of \k(S0 3 )p(0)(j{Di)\{C4))i, where 
\k{S0 3 )p{0)cf{Di)\{Ci))i are determined using the 
S0 3 D O D D 4 D C4 subduction. To do so, we use the 
relations given for k = 2 in Ref. 26 (p. 549): 



12211)3 = -^|2200) 4 + *i=|2222) 4 
|221-l) 3 = -i=|2200) 4 -z-^|2222) 4 
|2100)3 = (1 + i)-^=|2111) 4 

+ (l-i)-L|211-l) 4 --J=|2122) 4 



|2111>3 = (! + *)( 



1 



1 



2V2 2V6 
1 1 



)|21H>4 



+ ( - 1+t)( ivi + vf )|211 - 1>4 -vl |2122>4 

|2ii-i)3 = -(i + i) (^= + ^=)|2m) 4 

+ (1 - i)( ^-^i )|2il - 1) -> i ^ 

(41) 



If we now express the \k(S0 3 )p(0)a(D4)\(C4)}4 as 



In O, the d orbitals are expressed as: 



d xy = — (|22-2> - |22)) 
dy* = -^(|2-l> + |21)) 
dxz = ^(|2-l)-|21)) 



1% Z 2 _ r 2 



|20) 



d x 2_ y2 = —(\22) + \2-2)) 



(43) 



Combining Eqs 41, 42 and 43, we obtain the expres- 
sion of the d-functions, which are basis of the irreducible 
representations in C 3 , in function of the c?-orbitals in O: 



e +( e ) = |22H) 3 = -^d 3z 2_,.2 - -^d x 2_ y 2 



'.; ■. -■— r 2 "b ^d x 2_y2 



e_(e) = 1221-1)3 = ^ 
oi(t 2 ) = |2100)3 = 773 d ^ + 73 d ** - 73 d v* 
e+(i 2 ) = |2111>3 = 73^ + + ^ 

+ (i773 + 5)^ 
e+(i 2 ) = |211-1)3 = 77^ - (5775 + 



(44) 



In I?3, (e+(f 2 ), e_(t 2 )) is basis of the irreducible rep- 
resentation e(t 2 ). Similarly, (e + (e) , e+(e)) is a basis of 
the irreducible representation e(e). In D 3l a mixing is 
thus possible between the functions belonging to the two 
e irreducible representations, originating from the e and 
i 2 levels in O. 



F. Definition of the crystal-field parameters used 
in the LFM calculations 

In SO3 symmetry, the crystal-field Hamiltonian can 
be written as a combination of the q components of unit 
tensors with rank k. Each tensor is associated 
to the k irreducible representation of SO3. For the Cr 3+ 
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TABLE IV: Branching rules giving as irreducible represen- 
tation in D 3 



so 3 




O 


-> D 3 


-» c 3 




4 







-»■ 


-> 




4 




1 


-> 


-> 


j^-4100 


2 




1 


-> 








G. Relations between the crystal-field parameters 
used in LFM calculations and those derived from 
optical absorption spectroscopy 

As mentioned in Ref. 47, two equivalent parameter 
sets are available in optical absorption spectrocopy to 
describe the crystal-field in trigonal symmetry: (D a , D T , 
D q ) and (V, v\ D' ). In the following, we make connexion 
between the two sets. 



ion, k = 0, 2 or 4, the term k = contributing only to 
the average energy of the configuration. Again, we will 
not make use of the parity (± or g/u). To study the 
Cr 3+ ion in trigonal symmetry D 3 , we need to consider 
the subduction S0 3 dOd D 3 D C3. 

The crystal-field Hamiltonian is expressed as: 

ff cc = V" X HS0 3 )p(0)a(D 3 )\(C 3 ) u k(S0 3 )p(0)a(D 3 )\(C 3 )^ 

The unit tensor uKso 3 ) p {o)^d 3 )\(c 3 ) ig re i atec i to tne \ 
irreducible representation of C 3 subgroup, coming from 
the k irreducible representation of S0 3 , that becomes p in 
0, a in D 3 and A in C 3 . The terms x k( - s ° 3 ^°^ D3 ^ c ^ 
are the crystal-field parameters used in the LFM code. 
Their definition, in function of (D a , D T , D q ) or (y, v\ 
Dq), are given in Appendix G. 

The branching rules which give as irreducible represen- 
tation in D 3 are summarized in Tab. IV. This implies 
that : 



H rr = X 



4000 jj 4000 



+ X 4100 U 4100 + X 



2100^2100 



(45) 



1. (D a , D T , D q ) parameter set 



H cc is defined as H cc = H cub + H trig , where H cub and 
Htrig are the Hamiltonian for the cubic and the trigonal 
contributions to the crystal-field, respectively. According 
to Ref. 47 (Eq. 3.88), we have: 



(e±(t 2 )\Htrig\e±(t 2 )) =D a 



(ai(t 2 )\H trtg \ ai (t 2 )) = -2D a - 6D T , 
7 

(e±(e)\H trig \e±(e)) = -D T , 
(e±(t 2 )\H trig \e ± (e)) = ~^D a - 5D T ). 



If the cubic term H cu b is added, we have: 



(e±(* 2 )|ircc|e±(t2)) 



-AD q + D a + -D T , 



Using the expression of the d-eigenstates in C 3 given in 
Appendix E and the Wigner-Eckhart theorem (Eq. 4.2.4 
of Ref. 26), we can now calculate the matrix elements 

(k 1 (S0 3 )p 1 (0)a 1 (D 3 )X 1 (C 3 )\H cc \k 2 (S0 3 )p 2 (0)a 2 (D 3 )X 2 (C 3 )) 



(ai(t 2 )\H cc \ ai (t 2 )) = -AD q - 2D a - 6D T 
(e±(e)|ff cc |e±(e)) =6D q + -D T , 



(e±(t 2 )\H cc \e±(e)) 



V2 



(3D a - 5D T ). 



We obtain the following equations: 

V2 



{e±(t 2 )\H cc \e±(t 2 )} =- 



3V15 
1 



X 



4000 



V2 



3V21 



X 



4100 



Combining Eqs 46 and 48, we obtain: 
10 



2100 



(0l(t 2 )|Hcc|0l(t2)) 

(e±(e)|tf cc |e±(e)) 
(e±(i 2 )| J ff cc |e±(e)) 



70 

V2 



3V15 



X 



4000 



n V2 



^4000 _ 



^4100 
j^2100 



(18D, + 7D T ), 

Dr, 



V3 
V70D CT . 



3V21 



+ 



-2100 



(47) 



(48) 



(49) 



30 
1 



70 

j^4000 



2V2T 



X 



41000 



35 



X 



2100 



(46) 



X 4000 contains only the cubic part of the crystal field, 
although D T appears in its expression. This means that 
the trigonal distortion, via D Tl contributes to the e — t 2 
splitting of the d-orbitals. On the contrary, X 4100 and 
A" 2100 are entirely due to the trigonal distortion. Hence, 
when we investigated the effect of the trigonal distor- 
tion on the XANES spectra in the LFM calculations, the 
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value of X 4100 and X 2100 were set to zero, while the value 
of X 4000 was fixed to the value used in D 3c [ symmetry. 
Things can be simplified by defining D' ql so that it con- 
tains also the contribution of D T to the cubic field. This 
leads to the definition of two other distortion parameters, 
v and v' . 



2. (v, v' , D' q ) parameter set 
The crystal-field Hamiltoninan H cc is now defined as 



■HI 



trig 



where H' cub contains the contribution 



of D T . According to Ref. 47 (Eq. 3.90), we have: 



Combining Eqs 46 and 51, we have: 



X 



4000 



X 



4100 



6V30D' q , 
2V3 



V7 



x 2ioo = 4V35 1 
7 2V2 



3. Relations between the two parameter sets 



(52) 



I Combining Eqs 49 and 52, we obtain the relations given 

{e±{t 2 )\H' trig \e±{t 2 )) = --v, i n R e f. 47 (Eq. 3.91): 

2 

(ai(i2)|i? tr . is |ai(^)) = gi/, ? 

<e ± (i 2 )|i^|e ± (e)) = u' . D' q = D q + -D T , 

( 5 °) v = -3D, - f Dr , 

If we add the cubic term H' cub , which is here defined so f _ y/2 

that the center of gravity is maintained for the trigonally v = v + 5-^— D T . 
split t 2g orbitals, the following equations are obtained: 



1 



q 3 



(e±(t 2 )\H cc \e±(t 2 )) = W g 

2 

-; 

3 

(e±(i 2 )|flcc|e±(e)> =v' . 



( ai (t 2 \H cc \ ai {t 2 )) =-4D' q + -v, 



(51) 



(53) 



